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SUMMARY 

In  HP  turbines,  predictions  of  the  heat  transfer  to  the  blade  and  endwalls  is 
panicularly  important  for  an  accurate  assessment  of  turbine  component  life.  On  the 
endwalls,  there  are  often  complex  3D  (secondary)  flows  present  which  make  predictions  of 
heat  transfer  particularly  difflcult 

A  detailed  investigation  of  this  area  has  been  carried  out  on  a  fully  annular  cascade  of 
highly  3D  nozzle  guide  vanes.  Measurements  were  made  on  the  vane  and  endwalls  to 
determine  heat  transfer  and  aerodynamic  characteristics.  Testing  was  conducted  in  a  short 
duration  Isentropic  Light  Piston  test  facility,  at  engine  representative  Reynolds  number, 
Mach  number  and  gas-to-wall  temperature  ratio.  Interpreted  test  data  are  compared  with 
computations  obtained  at  test  conditions. 
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ABSTRACT 

In  HP  turbines,  predictions  of  the  heat  transfer  to 
the  blade  and  endwalls  is  pardculariy  important  fur  an 
accurate  assessment  of  turbine  component  life.  On 
the  endwalls,  there  are  often  complex  3D  (secondary) 
flows  present  which  make  predictions  of  heat  transfer 
particularly  difficult. 

A  detailed  investigadon  of  this  area  has  been  carried 
out  on  a  fully  annular  cascade  of  highly  3D  nozzle 
guide  vanes.  Measurements  were  made  on  tne  vane 
and  endwalls  to  determine  heal  transfer  and 
aerodynamic  characteristics.  Testing  was  conducted 
in  a  short  duration  Isentropic  Light  Piston  test 
facility,  at  engine  representative  Reynolds  number, 
Mach  number  and  gas-to-wall  temperature  ratio. 
Interpreted  test  data  are  compared  with  computations 
obtained  at  test  conditions. 

I,  INTRODUCTION 

In  order  to  increase  thrust*to-weight  ratio  and 
achieve  maximum  cycle  efficiencies  with  gas  turbine 
engines  it  is  necessary  to  raise  the  cycle  temperatures 
to  the  maximum,  within  the  constraints  of  structural 
integrity.  Thus  the  need  to  understand  in  detail  and 
predict  accurately  the  heat  transfer  distributions  for 
high  pressure  turbines  becomes  an  important  factor. 
The  presence  of  complex  highly  thrw-dimensional 
secondary  flows  within  the  turbine  passage  makes  the 
turbine  designer’s  task  very  difficult. 

Several  different  investigations  of  nozzle  guide  vane 
(NGV)  aerofoil  and  endwall  heat  transfer  behaviour 
have  been  repotted  in  the  literature.  The  work  of 
York  el  al  (1983),  Gladden  et  al,  (1988)  and  Boyle  et 
al  (1989),  concentrated  on  two-dimensional  cascade 
measurements,  and  although  general  features  of  the 
flow  were  modelled  radial  pressure  gradients  were 
not.  Granziani  et  al  (1979),  and  Caugler  and  Russell 
(1983),  correlated  detailed  heat  transfer 
measuRments  with  secondary  flow  patterns.  This 
work  was  not  represenutive  of  engine  operating 
conditions  in  temts  of  Reynolds  number,  Mach  number 
and  gas-to-wall  temperatuie  ratio.  Mote  recent  work 
by  Harasgama  et  al,  (1990),  and  Harvey  et  al,  (1990), 
have  shown  detailed  information  on  the  pattern  of  heat 
transfer  within  a  NGV  passage.  However,  these  NGV 
designs  incorporated  relatively  few  three-dimensional 
features. 


A  highly  three-dimensional  nozzle  guide  vane  has 
been  designed  and  tested  at  DRA  Pyestock.  The 
prime  objective  was  to  produce  a  challenging 
three-dimensional  design  which  would  act  as  a 
stimulus  for  aerodynamic  and  thermodynamic 
research.  Fully  three-dimensional  flow  calculation 
rmthods  were  employed  in  the  design  of  the 
variable-lean  nozzle  guide  vane.  Figure  1  shows  the 
details  of  the  vane  and  the  computational  grid  used. 

The  purpose  of  this  paper  is  to  describe  heat  flux  and 
pressure  measurements  obtained  from  a  set  of  NGV’s 
using  the  Isentropic  Light  Piston  Facility  (ILPF)  a 
short  duration  heat  transfer  tunnel.  The  NGV 
annular  nng  was  tested  over  a  range  of  Reynolds 
number,  .Mach  number  and  gas-to-wall  temperature 
ratios,  representative  of  engine  conditions.  The 
detailed  measurements  are  compared  with  theoretical 
flow  and  heat  transfer  predictions. 


Figure  I.  HTDU  4X  vane  and  compulalionni  grid 


2  DESIGN  DETAILS 
2.1  Duty 

The  aerodynamic  design  of  the  high  pressure  turbine 
was  targeted  for  a  future  engine  of  advanced  military 
duty.  The  turbine  has  been  designated  the  High 
Temperature  Demonstrator  Unit  (HTDU)  4X.  It  has 
the  following  overall  design  parameters. 

Stator  outlet  temperatuie  (SOT)  2066. 1 K 

Specific  work  capacity  (QpAT/T)  179.6  J/Kg  K 

Stage  loading  factor  AH/U^  2.002 

Stage  pressure  ratio  2.416 

Flow  function  VaAl*  0.625 

NGV  flow  function  wi/T/P  1.06  x  10  ’  KgVK  mVijj^j 
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converted  to  heat  transfer  rate  by  an  electrical 
analogue  circuit  which  solves  the  one-dimensional 
transient  heat  conducdon  equation  (Oldfield  et  al, 
1984).  The  output  signal  is  then  digitised  and 
recorded  on  a  nuni  computer.  The  present  tests 
included  thin  film  gauges  at  10%,  20%.  40%,  60%,  80% 
and  90%  spanwise  locations  with  12  on  the  suction  side 
and  10  on  the  pressure  side  at  each  spanwise  location. 
There  were  82  gauges  on  the  inner  w^  and  76  gauges 
on  the  outer  wall.  With  this  level  of  instrumentation 
it  was  possible  to  obtain  an  accurate  distribution  of 
vane  surface  and  endwall  heat  transfer.  Accuracy  of 
the  present  heat  transfer  data  is  within  ±5%.  Figure  3 
shows  vanes  instrumented  with  thin  films. 


Figure  3.  Vanes  instrumented  with  thin  films 


The  distribution  of  static  pressure  on  the  NGVs  was 
measured  using  surface  tappings,  and  the  inlet  total 
pressure  was  detennined  using  3  inlet  probes, 
positioned  at  120°  intervals  circumferentially.  The 
pressure  signals  were  recorded  using  a  fast  acting 
Scanivalve’  ZOC  (Zero,  Operate,  (Mbrate)  system. 
The  NGVs  were  fitted  with  tappings  at  5%.  10%.  20%, 
40%,  60%.  80%.  90%  and  95%  span  with  14  on  the 
suction  side  and  11  on  the  pressure  side,  at  each 
spanwise  location.  The  inner  platform  was 
instrumented  with  83  pressure  tappings  and  the  outer 
platform  contained  77  pressure  tappings.  The  overall 
accuracy  for  the  aerodynamic  results  is  such  that  the 
isentropic  Mach  numbers  calculated  from  the  measured 
static  pressures  are  accurate  to  within  ±  0.1%  for 
transonic  flow  conditions.  At  low  flow  velocities  the 
error  is  significantly  larger,  representing  about  10% 
at  0.1  Mach  number  due  to  the  smaller  difference 
between  the  toul  and  static  pressures.  All  tests  were 
performed  with  a  turbulence  grid  at  4.5  axial  chords 
upstream  of  the  NGVs,  giving  an  inlet  turbulence 
level  of  6.5%.  Table  1  shows  details  of  the  NGV 
operating  conditions. 


4  AERODYNAMIC  RESULTS 
4.1  Surface  Aerodynamics 


little  influenced  by  changes  in  exit  Mach  number. 
The  suction  side  results  are  similar  for  all  three 
conditions  during  the  initial  acceleration  from  the 
leading  edge  to  approximately  20%  axial  chord.  The 
below  design  and  design  Mach  number  condition  show 
diffusion  once  the  initial  acceleration  peak  has  been 
reached,  whereas  the  above  design  condition 
continues  to  accelerate  further  before  diffusion  takes 
place  at  approximately  85%  axial  chord. 

Measured  and  predicted  Mach  number  distributions  at 
60%  span  at  the  design  condition  are  shown  in  Figure 
5.  As  can  be  seen  the  Dawes  code  (Dawes,  1986) 
does  not  predict  the  suction  side  Qiffcsion  very  well. 
Possible  reasons  for  this  could  be  inadequate 
modelling  of  the  trailing  edge  geometry,  differe.nces 
between  the  modelled  and  true  downstream  annulus 
duct  and  the  lacx  of  grid  resolution.  The  grid  used 
for  this  prediction  is  74  axial,  25  radial  and  25 
tangential  points  on  a  sheared  H  grid.  The 
discrepancies  near  the  leading  edge  on  the  pressure 
side  are  attributed  to  experimental  error.  The  total 
to  static  pressure  differences  are  small  in  this  region, 
thus  the  error  in  isentropic  Mach  number  becomes 
more  significant.  The  remainder  of  the  pressure 
surface  is  very  \ve:l  modelled. 
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Figure  4.  60%  span  Mach  number  variation 


Mach  number  distributions  at  60%  span  for  three  Mach 
number  conditions  are  shown  in  Figure  4.  The 
measured  Reynolds  number  has  been  kept  ccnsiant  for 
these  results.  The  pressure  side  aer^ynamics  are 
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Measured  and  predicted  suction  side  isentropic  .Mach 
number  distributions  for  three  exit  .Vlach  number 
conditions  are  shown  in  Figure  6.  Mach  numbers  at 
the  vane  root  are  higher  than  the  up.  as  required  by 
radial  equilibrium.  The  distributions  indicate  a 
significant  change  in  aerodynamics  as  the  exit  Mach 
number  is  increased  from  the  below  design  condition 
to  design  and  above  design.  The  predictions  show 
similar  features  to  the  experimental  results.  The 
initial  rapid  acceleration  is  predicted  well  in  all  three 
cases,  with  some  differences  towards  the  tip. 


may  have  introduced  small  local  differences  in  these 
areas. 


Measured 

Predicted 
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Figure  5.  Vane  surface  Mach  number 
distribution  at  60%  span. 

However,  this  region  has  highly  3D  geometry  and  is 
of  particular  difficulty  for  instrumenting  static 
tappings.  It  is  possible  that  some  tappings  may  not 
be  normal  to  the  local  surface.  The  vanes  are 
manufactured  by  the  use  of  a  computer  numencally 
controlled  machine  and  are  then  hand  finished  to 
achieve  a  smooth  surface  particularly  at  the  leading 
edge,  at  the  junction  between  endwalls  and  the  vane 
and  at  the  trailing  edge.  This  final  finishing  process 


Measu 


Discrepancies  towards  the  trailing  edge  of  the  vane 
vary  depending  on  the  Mach  number  condition.  For 
the  below  design  case  the  prediction  compares  very 
well  with  a  few  small  differences  near  the  trailing 
edge  berween  40%  and  60%  span.  As  previously 
discussed,  the  design  and  above  design  conditions 
show  much  less  agreement  with  the  prediction  once 
Che  inidai  acceleration  peak  is  reached.  The 
experimental  results  show  mainly  diffusion  but  the 
prediction  indicates  a  further  acceleration. 

The  measured  pressure  side  iseniropk  Mach  number 
for  the  design  condition  is  shown  in  Figure  7  together 
with  a  predicnon.  The  prediction  compares  very  well 
with  the  computed  results.  This  was  also  true  of  the 
below  and  above  design  Mach  number  condition. 

Hub  and  casing  .Mach  number  distributions  for 
measured  a.nd  predicted  results  are  shown  in  Figure 
8a  and  8b.  The  hub  Mach  numbers  exceed  those  on 
the  casing  due  to  the  radial  equilibrium  criteria. 
When  operating  at  the  high  Mach  number  condition  the 
hub  platform  shows  a  trailing  edge  shock  impinging  on 
the  suction  surface.  The  compa."isQn  between 
measured  a.nd  computed  results  shows  very  good 
agreement,  for  both  hub  and  casing  platforms.  These 
results  obviously  validate  the  3D  Navier-Siokes  code 
as  a  f.ist  and  relatively  accurate  rreihod  for  the 
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Individual  experimental  sweeps,  after  processing, 
have  been  combined  to  form  the  complete  area 
traverse.  The  traverse  was  carried  out  at 
approximately  5%  axial  chord  downstream  of  the 
trai  ing  edge.  The  measured  profile  wakes  can  be 
clearly  seen,  at  mid-height  the  wake  is  considerably 
wider  than  that  near  the  endwalls.  The  wake  is  wider 
due  to  the  trailing  edge  stacking;  the  probe  tip  is 
thus  further  from  the  trailing  edge  at  mid-height  than 
near  the  root  and  tip.  The  measured  data  also  show 
areas  of  significant  total  pressure  loss,  these 
coiiespond  to  the  horse-shoe  vortex  as  seen  in  Figure 
12a.  The  Navier-Stokes  prediction  also  sliows  wider 
wakes  at  mid-height  than  netu  the  endwalls  and 
compares  very  well,  although  the  tueas  of  much  larger 
loss  created  by  the  horse-shoe  vortex  are  not 
predicted  at  all. 

5  HEAT  TRANSFER  RESULTS 

Nusselt  number  distributions  on  the  vane  aerofoil 
surface  at  60%  span  for  3  running  conditions  are 
shown  in  Figure  10.  As  can  be  seen  the  .Nusselt 
number  increases  as  Reynolds  number  is  increased, 
on  both  suction  and  pressure  side.  The  variation  of 
heat  transfer  with  Mach  number  is  also  shown  on 
Figure  10.  In  this  case  Nusselt  number  decreases 
with  increasing  Mach  number.  This  is  due  to  the 
shear  stress  in  the  boundary  layer  decreasing  as  the 
flow  velocity  increases.  These  results  are  in 
agreement  with  the  work  of  previous  authors  (Jones 
etal,  1989). 

Contours  of  Nusselt  number  overlaid  onto  a  3D  view  of 
the  NGV  suction  side  are  presented  in  Figure  11.  it 
can  be  seen  that  there  is  a  significant  effect  of  the 
secondary  flow  on  the  heat  transfer  distributions. 
The  effect  is  seen  to  be  more  pronounced  at  the  tip 
than  the  hub.  The  mid-span  Nusselt  numbers 


towards  the  trailing  edge  are  consistently  higher  than 
near  the  endwalls,  for  all  operating  conditions.  This 
is  due  to  endwall  boundary  layer  migration  onto  the 
vane  suction  side  resulting  from  secondary  flow  in  the 
vane  passage.  Flow  visualisation  performed  on  the 
endwall  and  vane  surface  in  the  ILPC  (Figure  12a  and 
I2b)  indicate  that  the  secondary  flow  is  stronger  near 
the  casing  than  at  the  hub.  This  effect  is  well 
predicted  b>  liie  Dawes  (1986)  How  solver  as  can  be 
seen  in  Figure  13.  The  hub  secondary  How  appears 
on  the  suction  surface  later  than  the  casing  and  docs 
not  migrate  along  the  span  to  the  same  extent.  The 
prediction  shows  this  effect  but  to  somewhat  lesser 
extent  than  indicated  by  the  flow  visualisadon. 

Increasing  Reynolds  number  causes  the  Nusselt 
number  over  the  entire  vane  surface  to  increase, 
whilst  on  the  suction  surt'ace  ihe  reverse  occurs  for 
increasing  Mach  number.  The  pressure  side  Nusselt 
number  distributions  did  not  vary  significantly  with 
the  Mach  number  changes  and  there  is  very  little 
span  wise  variation,  see  Figure  14. 

Endwall  heat  transfer  is  affected  by  the  inlet 
boundary  layer,  vane  aerodynamic  loading  and  the 
geometry  of  the  vane  passage.  The  tests  have  shown 
that  there  are  some  di^erences  in  the  inner  and  outer 
wall  Nusselt  numbers  both  in  magnitude  and 
distribution.  Figure  13  shows  hub  Nusselt  numbers 
for  all  conditions.  As  the  Reynolds  number  is 
increased  the  region  of  highest  heat  transfer  becomes 
concentrated  at  the  pressure  side  trailing  edge. 

The  inlet  endwall  boundary  layer  stagnates  at  the 
vane  leading  edge  and  forms  a  horseshoe  vortex.  The 
two  legs  of  this  vortex,  on  the  pressure  and  suction 
sides,  proceed  into  the  two  adjacent  passages.  The 
pressure  side  leg  is  swept  across  the  endwall  by  die 


Figure  8b.  Measured  and  predicted  casing  Mach  number  distribution 


passage  secondary  vortex  to  the  suction  side  and  then 
up  ihe  vane  suction  surface.  Tliis  causes  the 
boundary  layer  to  be  stripped  off  the  endwall  and 
convected  to  the  suction  side  of  the  passage  and  then 
up  the  suction  surface  of  the  vane.  A  new  thin 
highly  skewed  boundary  layer  thus  forms  behind  the 


separation  line,  which  results  in  high  heat  transfer  in 
this  region  (Ghana,  1992).  Low  values  of  heat 
transfer  are  evident  close  to  the  suction  side  of  the 
passage  after  approximately  50%  axial  chord.  This  is 
generated  by  the  low  energy  fluid  migrating  across 
from  the  pressure  side  trailing  edge  region.  Outer 
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Figure  9.  Measured  and  predicted  total  pressure  contours  105%  axial  chord 


Figure  12.  Suction  side  flow  visualisation 
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Figure  14.  Pressure  side  Nusselt  number 
distribution 

An  attempt  has  been  made  to  compute  the  Nusselt 
number  distnbution  to  the  vane  at  &07o  span.  Figure 
17  shows  measured  and  predicted  aerofoil  Nusselt 
numbers.  The  predictions  were  performed  using 
versions  of  the  Stancool  and  Texstan  boundary  layer 
codes  (Crawford  and  Kays,  1976).  The  codes  were 
run  using  the  predicted  pressure  distribution  from  3D 
Dawes.  Two-dimensional  calculations  were  performed 
along  grid  lines  on  the  vane  surfaces;  no  allowance 


for  surt’acc  curvature  was  included.  The  gr.d  lines  in 
this  region  follow  the  stream  lines  very  closely.  The 
inlet  turbulence  was  set  at  6.59b.  The  Prandtl  mixing 
length  mocei  was  used  with  Van  Driest  damping  in  the 
near  wall  region. 

The  Texsian  computation  was  stance  close  to  the 
leading  edge  from  a  specified  boundary  layer 
dimension,  in  the  case  of  Stancool  an  initial  distance 
from  the  stagnanon  point  was  specified.  It  was  found 
that  the  subsequent  turbulent  boundary  layer 
development  was  only  weakly  dependent  on  the  initial 
laminar  region.  Transition  location  was  determined  by 
the  correlation  of  Abu-Ghannam  and  Shaw  (1980)  in 
the  case  of  Stancool  computations  and  a  transition 
Reynolds  number  of  200  was  specified  for  Texstan.  It 
will  be  seen  that  l.'■ansition  is  quite  well  predicted 
using  both  approaches.  The  heat  transfer  levels  are 
in  error  by  around  309b  in  me  worst  case,  for  the 
pressure  side,  out  the  errors  are  more  typically  10%. 
The  pressure  surface  results  from  Texstan  show  good 
agreement.  The  discrepancies  towards  the  sailing 
edge  are  probabiy  due  to  irregular. ties  in  the  trailing 
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Figure  15.  Hub  Nusselt  number  distributions 
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edge  predicted  pressure  distribution  from  Dawes,  as 
discussed  in  Section  4.  On  the  suction  surface 
Texstan  performed  less  well  with  a  large  over-shoot 
but  then  approached  the  experimental  results. 
Stancool  predicted  the  suction  side  Nusselt  number 
distribution  reasonably  well  but  substantially  under 
predicted  on  the  pressure  side.  Stancool  .xsulis  with 
the  transition  location  determined  by  a  Reynolds 
number  of  2(X)  predicted  slightly  later  transition  and 
gave  slightly  better  results. 


remain  tc  be  explored.  These  include  smoothing  of 
the  imposed  free  stream  velocity  distribution, 
curvature  specification  and  introduction  of  turbulent 
iranspon  equations. 

CONCLUSIONS 

A  highly  j-dimensional  nozzie  guide  vane  has  been 
successfully  designed  ,md  tested  at  engine 
xprcscntative  conditions. 


The  Texstan  code  used  was  an  earlier  release  and  did 
not  include  many  of  the  features  of  the  more  mature 
Stancool  code.  Many  refinements  to  the  predictions 


Three  sets  of  aerodynamic  measurements  have  been 
taken  at  different  Mach  number  conditions.  The 
results  indicate  the  vane  is  generally  operaung  as 
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Figure  16.  Casing  Nusscit  number  distributions 
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designed  with  some  unexpected  behaviour  towards  the 
suction  side  trailing  edge.  This  warrants  further 
investigation.  The  flow  visualisation  of  endwall  cross 
flow  indicated  the  presence  of  secondary  tlows  and 
ties  in  with  measured  heat  transfer  patterns.  Regions 
of  high  heat  transfer  found  in  other  designs 
previously  (Harasgama,  1990)  in  the  endwail  ;nid 
passage  region  have  been  reduced  in  lU'ea  and 
confined  to  the  trailing  edge  region.  Variation  of  the 
Nusselt  number  with  changes  in  Reynolds  number  are 
consistent  with  classical  flat  plate  analogy.  The 
effect  of  Mach  number  vanation  on  Nusselt  number 
confirms  results  reponed  by  other  authors. 
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All  the  measurements  and  predictions  indicate  that  the 
secondary  flows  present  are  strongest  near  the 
casing.  The  downstream  area  traversing  technique 
has  produced  good  results,  and  ihese  compare  well 
with  those  predicted  using  the  3D  Dawes  Navier 
Stokes  solver.  However  the  computations  made  so  f.u- 
show  no  evidence  of  the  endwall  horse-shoe  voruces. 

Boundary  layer  heat  transfer  calculations  at  60%  span 
on  the  vane  surfaces  are  promisf-.g.  Transition  was 
well  predicted  using  either  a  Reynolds  number  of  200 
or  the  correlation  of  Abu-Ghannam  and  Shaw. 
Predicted  heat  transfer  levels  .ue  in  error  by 
typically  10%  with  a  worse  case  error  of  30%  Further 
work  is  obviously  required  to  resolve  inconsistencies 
and  refine  the  results. 


-  aa--  -  - _ _ _ _ 
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Figure  17.  Predicted  and  measured  aerofoil 
Nusselt  number  at  60%  span 
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NOMENCLATURE 

Cp  -  Specific  heal  at  constant  pressure 

AH  -  Enthalpy  change 

Mn  ■  Isentropic  Mach  number 

Re  -  Reynolds  number  based  on  true  chord 

P  -  Total  pressure 

T  -  Total  temperature 

AT  -  Temperature  change 

Va  -  Axial  Velocity 

U  -  Mean  Blade  Speed 
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